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RESUMEN 

El resumen sera traducido al espanol por los editores. In dense stellar 
systems the frequent dynamical interactions between stars play a crucial role in 
the formation and evolution of compact binaries. We study these processes using a 
novel approach combining a state-of-the-art binary population synthesis code with a 
simple treatment of dynamical interactions in dense star cluster cores. Here we focus 
on the dynamical and evolutionary processes leading to the formation of compact 
binaries containing white dwarfs in dense globular clusters. We demonstrate that 
dynamics can increase by factors 2 — 100 the production rates of interesting 
binaries such as cataclysmic variables, "nonflickerers" (He white dwarfs with a 
heavier dark companion), merging white dwarf binaries with total masses above 
the Chandrasekhar limit, and white dwarf binaries emitting gravitational waves in 
the LISA band. 

ABSTRACT 

In dense stellar systems the frequent dynamical interactions between stars 
play a crucial role in the formation and evolution of compact binaries. We study 
these processes using a novel approach combining a state-of-the-art binary popu- 
lation synthesis code with a simple treatment of dynamical interactions in dense 
star cluster cores. Here we focus on the dynamical and evolutionary processes lead- 
ing to the formation of compact binaries containing white dwarfs in dense globular 
clusters. We demonstrate that dynamics can increase by factors ^ 2 — 100 the pro- 
duction rates of interesting binaries such as cataclysmic variables, "nonflickerers" 
(He white dwarfs with a heavier dark companion), merging white dwarf binaries 
with total masses above the Chandrasekhar limit, and white dwarf binaries emitting 
gravitational waves in the LISA band. 
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1. INTRODUCTION 

From the earliest observations of X-ray binaries 
in globular clusters it has always been clear that they 
must be very efficient fact ories for the production 
of compact binary systems llClarklh 97,^1^ . The over- 
abundance of compact binaries in clusters, as com- 
pared to the field, must be a result of close stellar 
encounters. The key processes that affect the binary 
population in dense cluster environments include the 
destruction of wide binaries, ha rdening of clos e bina- 
ries (following "Heggie's law" (Hcggie l ll975j) V and 
exchange interactions, through which low-mass com- 
panions tend to be replaced by a more massive par- 
ticipant in the encounter. As a result of these pro- 



cesses, in the dense cores of globular clusters, bi- 
naries are strongly depleted and their period distri- 
bution is very different from that of a field popu- 
lation (Ivanova et al. 2004). These processes also 
lead to an interesting and complex interplay between 
dynamics and binary evolution. For example, ex- 
change interactions involving compact objects often 
produce systems that will evolve through a common- 
envelope (CE) phase and form very short-period bi- 
naries j_wlnch_ar£jnuchj^^ popula- 
tions ijRasio. Pfahl. fc RaDPapOT3l2000() . 

There are two possible approaches to the 
study of binary evolution and dynamics in glob- 
ular clusters. One can start from A'^-body sim- 
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ulations and introduce various simplified treat- 
ments of binary star evolution. This has 
been the traditional approach f or many years 
(for recent examples see Ishara fc Hurlevl l2002t 
IPortegies Zwart. McMillan. Hut, fc Makind l200l|) 
Alternatively, one can start from a binary popula- 
tion synthesis code and add a treatment of dynam- 
ical interactions. This approach was pioneered by 
Portegies Zwart et al. (1997) and has been adopted 
in our recent work. It has the great advantage that it 
is computationally much less expensive than TV-body 
simulations, so that more exploration of the (enor- 
mous) parameter space is possible, and more realistic 
simulations, using sufficiently large numbers of stars 
and binaries, are possible today. In contrast, even 
when using special-purpose GRAPE computers, A''- 
body simulations are still limited to smaller systems 
like open clusters with limited coverage of parameter 



space and with unrealisticalb 


r small numbers of bina- 


ries (seellvanova et al. 1 


I2nn4l 


IWilkinson et al.ll200.'^. 



In our code we combined StarTrack, a 
state-of-the-art binary population synthesis code 
l)Belczvnski et al.l l2002l) and FewBody, a small- 
A''-body integrator that we use to compute 3- 
body and 4-body interaction s ((Fregeau et al. lEf)o3 
iFregeau fc RaDPaportll2004(l . Currently we adopt 
a simple two-zone model, in which the clus- 
ter is partitioned into an inner core and an 
outer halo, with all interactions assumed to take 
place in the core. This background cluster 
model remains unchanged throug hout the evolution 
ifHut. McMillan, fc Romanil99a) . In particular, the 
core density is assumed constant. However, our ul- 
timate aim is to incorporate full dy namical Monte 
Carlo models ( Fregeau et al. lE'00.1^ . In a typical 
simulation we start with N ~ 10^ stars, with be- 
tween 50% and 100% binaries. This high primor- 
dial binary fraction (much higher than assumed in 
all previous studies) is needed in order to match the 
observ ed binary fractions i n globular cluster cores 
today Ijlvanova et al. 



2. COMPACT BINARIES AND MERGERS 
Consider the evolution of a typical dense cluster 
with central velocity dispersion a = lOkms"^ and 
core density ric = lO^pc"^. Dynamical interactions 
lead to greatly enhanced numbers of compact bina- 
ries containing white dwarf (WD) and neutron star 
(NS) components, and, in particular, to much larger 
numbers of heavier compact binaries (Fig. 1). Here 
we focus in particular on the fate of WDs involved 
in CE events leading to compact binary formation or 
mergers of WDs driven by gravitational wave emis- 
sion. Our simulations confirm that the formation 
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Fig. 1. Mergers and dynamical formation of compact 
binaries containing white dwarfs (WD) or neutron stars 
(NS). The horizontal axis gives the cluster age when the 
event occurred and the vertical axis gives the total mass 
of the binary. The top panel shows a population evolved 
without dynamical interactions, while the bottom panel 
shows the same initial population evolved in a dense clus- 
ter with core density ric — lO^pc"^. The size of each 
symbol indicates the types of stars involved, while the 
shape indicates the type of event: binary merger (round 
dots; usually driven by gravitational radiation or spiral- 
in during a common envelope phase), compact binary 
formation following a common envelope phase (stars), 
or merger following a physical collision (triangles). It 
is clear that dynamics increases both the total number 
of events and the typical mass of mergers and compact 
binaries. 
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rate of compact binaries containing a Helium WD 
with a heavier companion via CE events is increased 
significantly in dense clusters. The brighter Helium 
WDs in these binaries could be detectable as "non- 
flickerers." These were o bserved for the first time 
in the core of NGC 6397 l)Cool et al. 199^) and are 
indeed thought to be double WD binaries containing 
a young Heliu m WD with an older and hea vier WD 
companion (H ansen. Kalogera. fc Rasiol |20031 . 

We can also examine the rate of double WD 
mergers, and, in particular, those for which the to- 
tal mass is > Mch — IAMq). These so-called 
supra-Chandrasekhar mergers could lead either to a 
Type la supernova, or to a "merger-induced" col- 
lapse of the remnant to form a NS (and perhaps 
a millisecond radio pulsar). It is possible that the 
NS in this case is formed without a significant kick 
(see the article by Podsiadlowski in this volume). 
An increased rate of Type la supernovae from star 
cluster dynamics would likely be redshift-dependent 
(as more stars are formed in starburst environments 
— which favor star cluster formation — at higher 
redshifts) and this has important potential conse- 
quences for their use in cosmology (for a review see 
iLeibundeutj l200l|) . Alternatively, merger-induced 
collapse could lead to the formation of neutron stars 
and millisecond pulsars in clusters, thereby alleviat- 
ing or perhaps eliminating the NS "retention prob- 
lem" (the ejection of most NS from the shallow clus- 
ter potentials if they are born with the natal kicks 
expected from asymmetric supernova explosions; see, 
e.g., Chen & Leonard 1993; Pfahl, Rappaport, & 
Podsiadlowski 2002) 

The enhanced production rate of double WD 
mergers in dense stellar clusters w as first discussed 
in detail bv lShara fc HurlevI l)2002() . They estimated 
that, for stars born in open clusters (which can be 
simulated directly using their A^-body code), the 
supra-Chandrasekhar WD merger rate can be in- 
creased by an order o f magnitude. The results of 
IShara &: HurlevI ll2002l) are based on A^-body simu- 
lations for a typical open cluster containing ~ lO'' 
stars with 10% primordial binary fraction, and with 
a = 2kms~^ and — lO^pc"'^. Our simula- 
tions are for much denser and massive star clusters 
with a primordial binary fraction of 100%. Not sur- 
prisingly, we find an even larger number of supra- 
Chandrasekhar WD mergers, with an enhancement 
factor (compared to a field population) closer to 
^ 100 for a typical dense globular cluster. In con- 
trast, the total number of double WD mergers (of 
all types) is typically increased by a factor of a few 
only. The majority of these mergers are driven ul- 
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Fig. 2. All LISA-detectable binaries during the last Gyr 
of evolution for a typical dense cluster model with core 
density ric — pc~'^. Each source is shown only once, 
at a randomly chosen moment in its evolution across the 
LISA band. Here P is the binary orbital period and 
Mchirp = {MiM2f^^/{Mi + M2f^^ is the chirp mass. 
The different symbols indicate binaries with different 
types of components. We again compare a field pop- 
ulation (top) to our cluster model including dynamical 
interactions (bottom). 
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timately by gravitational radiation, although a few 
come from physical collisions of WDs during hard 
binary encounters (Fregeau et al. 2004). 

3. LISA SOURCES 

Our results show that compact double WD bina- 
ries are mainly formed dynamically and have typi- 
cally experienced multiple hardening encounters be- 
fore merging. Prior to merger, these systems may be 
detectable as gravitational wave sources by LISA, 
when their orbital period becomes smaller than 
about 2000 s (Benacquista, Portegies Zwart & Rasio, 
2001). This limit on the orbital period is imposed by 
the background noise from Galactic binaries. At the 
same time, the positional accuracy of LISA is much 
greater for binaries with these shorter periods, so 
that the sources can then be associated with specific 
globular clusters in our Galaxy. 

Figure 2 shows all LISA-detectable sources (chirp 
masses and periods) that appeared during the last 
Gyr of our simulation for a typical dense cluster. 
For this model, where the total cluster mass today 
is about 2 x 10^ Mq, at least one LISA source is 
present at any given time. On average, there are 5 
LISA sources at any given moment, about twice the 
number predicted for a field population. In addition 
to this rise in the number of sources, we also note 
changes in their typical properties: in particular, the 
number of sources with chirp mass above 0.4 Mq is 
increased significantly in the cluster model. We also 
find that NS-WD binaries represent about 20% of 
all LISA sources. Therefore, the number of LISA- 
detectable NS-WD binaries in all Galactic globular 
clusters (Mtot ^ 10 ''■^Mq) could be as high as > 
100, while the total number of detectable WD- WD 
binaries could be > 500. However, since the cluster 
models we have used so far in our simulations are 
denser than average, these numbers should be taken 
as upper bounds. 

4. MASS TRANSFER SYSTEMS 

As they spiral in and evolve across the LISA 
band, WD- WD and NS-WD binaries will eventually 
come into contact around a period P ~ 100 s. At 
first the binary orbit shrinks as it looses angular mo- 
mentum to gravitational radiation. However, during 
stable mass transfer, the orbit will evolve towards a 
larger period. During this mass transfer phase the 



binary can also appear as an ultracompact X-ray bi- 
nary (NS-WD) or an AM CVn type cataclysmic vari- 
able (WD- WD). As expected, we find that the popu- 
lation of these mass transfer binaries is also increased 
significantly by dynamical interactions. For exam- 
ple, for our typical cluster model with ric = 10^ pc~^, 
we predict about 50 AM CVn binaries, roughly twice 
the number obtained without dynamics. 
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